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NOMENCLATURE 


(X  -  Lift  ctirve  slope: 

b  -  Blade  semichord  at  each  radial  station. 

K  -  »*Effective"  plunging  velocity  of  the  blade  section; 

includes  the.,  .component  of  the  free  stream  velocity,  , 
normal  to  the  local  blade  sections  (i.e,,  the  effects 
of  the  local  blade  slopes). 

b  -  Subscript  indicating  position  in  the  rotor  disk. 

H  -  Lift  per  unit  span. 

-  Number  of  radial  bound  vortex  segments  representing 
blades . 

N  ^  -  Nmber  of  equally-spaced  azimuth  positions  at  \diich 
the  airloads  are  computed;  must  be  multiple  of  the 
number  of  blades  in  the  rotor. 

?  -  The  points  of  .the  rotor  disk  where  the  airloads  are 

computed. 

r  -  The.  radial  coordinate  of  the  rotor. 

R  -  Blade  radius. 

V  -  Tangential  velocity  in  tip-path-plane :  V  » Ilv-  G»  0(7 

\4  -  Free  stream  velocity  due  to  rotor  translation. 

UT  -  The  induced  velocity  component  normal  to  the  tip-path- 
plane  at  the  three-quarter-chord  of  the  section. 

o(  -  Instantaneous  geometric  angle  of  attack  with  respect 
to  the  tip-path-plane. 

/\ 

o<  -  Stalling  angle  of  attack  for  the  blade  sections. 
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0<«  -  Instantaneous  effective  angle  of  attack  with  respect 
to  the  tip-path-plane;  =  +  ur/v. 

o(j.  -  Longitudinal  tip-path-plane  inclination  with  respect 
to  the  free  stream  velocity,  Vf  , 

T  -  The  unknown  strength  of  the  bound  vortex  segment 
representing  the  blade  section, 

M  -  Advance  ratio;  Cos . 

9  -  Air  density. 

^  -  Azimuthal  position, 

XI  -  Rotational  speed  of  the  rotor. 


SUMMARY 


This  investigation  is  an  initial  effort  to  obtain  a  solution  to  the 
airloads  problem  for  a  rotating  wing  in  steacty--state  translational  flight. 

A  sufficient  amount  of  the  wake  detail  was  retained  to  enable  accurate 
computation  of  the  wake-induced  velocities  at  the  blade.  The  method  of 
solution  developed  is  numerical,  utilizes  high-speed  digital  coirqjutation, 
and  is  relatively  simple  to  use. 

Each  blade  is  represented  by  a  segmented  lifting  line  and  the  shed 
and  trailing  vorticity  distributions  are  represented  by  a  continuous  mesh 
of  segmented  vortex  filaments  originating  at  the  instantaneous  position  of 
the  blade.  The  wake  of  each  blade  of  the  rotor  is  different;  these  wakes 
change  with  the  azimuth  position  of  the  blades  to  correspond  to  the  instantan¬ 
eous  wake  configTiration  for  each  azimuth  position. 

The  method  developed  was  used  to  compute  the  rotor  airload  distri¬ 
butions  for  the  Be].l  HU-IA,  the  Sikorsky  H-3U,  and  an  NASA  model  rotor  in 
forward  flight.  The  computed  airloads  were  harmonically  analyzed  and  the 
radial  distribution  of  each  harmonic  compared  with  the  measured  distributions. 
In  general,  encouraging  agreement  is  obtained  in  these  con^jarisons.  Effects 
of  changing  some  of  the  parameters  of  the  computation  were  investigated. 
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imODUCTION 


The  problem  of  predicting  the  aeroelastic  response  of  rotating  wings 
is  generally  recognized  to  be  very  difficult.  The  industry  has  had  to  rely 
on  expensive  empirical  methods  and  past  experience  to  design  and  develop 
these  aircraft  components. 

It  is  believed  that  the  investigation  that  is  being  reported  herein 
can  be  put  in  its  proper  context  by  relating  it  to  the  continuing  research 
program  on  helicopter  dynamics  that  has  been  carried  out  for  many  years  at 
Cornell  Aeronautical  Laboratoiy, 

One  of  the  earliest  direct  anticedents  (initiated  in  19h^)  of  the  current 
effort  is  the  -Variable  Stiffness  Blade  Program",  Refs,  k  -  21.  That  study 
was  initially  aimed  at  determining  if  rotor  blade  structural  stiffhess  was 
an  iTOortant  factor  in  the  blade  fatigue  life.  Early  in  this  program  it  was 
found  from  strain  gage  measurements  (made  in  flight)  that  higher  harmonic  order 
bending  moments  of  large  amplitude  existed  when  a  multiple  of  the  rotor  ipm 
was  near  the  frequencj^  of  the  second  or  third  natural  bending  mode.  Thereafter, 
emphasis  of  the  program  was  concentrated  on  investigation  and  evaluation  of 
the  resonance  effects.  It  was  postulated  that  the  forcing  functions  exciting 
the  bending  modes  of  the  blades  (and  the  fuselage  modes)  arose  from  aerodynamic 
forces  whose  existence  could  be  explained  only  on  the  basis  of  nonuniform  inflow 
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velocities.  Comparison  of  f;eneraliz;ed  forces  derived  from  measured  strain 
data  with  generalized  forces  comouted  from  downwash  velocities  based  on 
Mangier' s  analysis  (Ref.  22)  tended  to  confirm  this  hypothesis.  It  was  also 
noted  that  Mangier' s  analysis  was  not  sufficiently  refined  to  permit  adequate 
blade  stress  calculations. 

Thus,  The  "H-?  Variable  Stiffness  Blade  Program"  resulted  in  the 
identification  of  the  complex-induced  velocity  distribution  as  the  aerodynamic 
part  and  the  resonance  conditions  as  the  structural  part  of  the  rotor  fatigue 
problem,  A  sendemplrical  approach  based  on  the  use  of  generalized  forces 
obtained  from  strain  measurements  was  recommended  for  analyzing  rotors.  It 
is  believed  that  this  approach  is  still  used  as  a  basis  for  the  analysis  of 
rotor  stresses. 

Recognition  that  the  rotor  bending  stress  problem  was,  in  fact,  an 
aeroelastic  problem  led  to  a  more  careful  consideration  of  the  possibilities 
of  rotor  flutter  and,  conversely,  the  possibilities  for  introducing  couplings 
between  blade  modes  that  would  alleviate  stresses.  Several  studies  character¬ 
ized  by  anal3’tical  developments  and  corroborating  experiments  were  made  (Refs, 
23  -  35).  It  was  noted  that  fair  agreement  between  theoretical  and  experi¬ 
mental  flutter  results  was  obtained  for  most  hovering  configurations  but  that, 
under  certain  conditions,  a  multiplicity  of  flutter  modes  was  obtained  •vdiich 
were  not  predicted.  The  multiplicity  could  not  be  accounted  for  on  the  basis 
of  mass  and  elastic  characteristics  alone,  and  it  was  concluded  that  the  wake 
strength  and  position  associated  with  the  blade  flutter  motions  v;as  a  major 
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effect.  In  fact,  the  name  "wake  excited  flutter"  was  anplied  to  certain  of 
the  flutter  modes  that  had  frequencies  which  were  nearlj'-  integer  multioles 
of  rotor  speed.  Further  analj^tical  work  resulted  in  establishment  of  a  two- 
dimensional  unsteady  aerodjniatnic  theory  for  a  hovering  rotor  (Ref,  28)  that 
predicted  the  multiple  flutter  loops  and  further  verification  was  obtained 
from  additional  flutter  and  response  tests  (Refs.  30  and  32).  Forvfard  flight 
flutter  investigations  (Refs.  and  35)  were  also  initiated,  although  it 
was  realized  that  no  suitable  unsteady  aerodjTiamic  theory  v;as  available  for^ 
this  flight  regime. 

The  results  reported  herein  can  thus  be  viewed  as  a  logical  development 
in  a  history  that  started  more  than  fifteen  years  ago.  It  began  with  the 
proper  identification  of  the  blade  oscillatory  stress  problem  and  the  develop¬ 
ment  of  an  interim  semiempirical  solution.  Then  it  proceeded  through  a 
systematic  series  of  aeroelastic  analyses  and  validating  tests,  which  led 
logically  to  the  present  effort. 

The  aeroelastic  problem  can  be  thought  of  as  consisting  of  two  parts: 
the  dynamic  response,  and  the  airloads.  They  are  interdependent  and  related 
in  that  the  airloads  cause  the  dj-namic  response  and,  at  the  same  time,  depend 
on  it.  Much  of  the  difficulty  in  obtaining  an  adequate  solution  to  the  aero¬ 
elastic  problem  is  in  the  airloads  part  of  the  problem  (specifically  in 
obtaining  an  adequate  aerodynamic  representation  of  the  blades  and  wake). 


Any  accurate  method  of  computing  the  airloads  must  adequately  predict 


the  Hahe-induced  velocities  at  the  blades  because  the  airloads  are  strongly 
Influenced  by  these  velocities.  However,  because  the  vortical  wake  of  the 
rotating  wing  is  extremely  com;)lex  and  difficult  to  adequatel;'-  represent 
mathematically,  the  ^.ractical  solution  of  the  aeroelastic  problem  has  been 
delayed,  Harl3'  attempts  to  solve  this  mroblem  anal3’t,icall;,’'  were  based  on 
relativel;.''  drastic  simnlifications  of  the  wake  of  the  rotor  to  make  them 
comnutationall^/  feasible.  The  modem  high-sneed  commuting  machines  of  todaj'" 
have  made  it  mossiMe  to  acco-mt  for  much,  more  of  the  detail  of  the  wake  than 
has  been  nossible  in  the  mast  and  thus  nermit  an  adequate  aerodynamic  remre- 
sentation  of  the  blades  and  wake  to  be  formulated  which  vrill  enable  accurate 
corrriutation  of  the  rotor  airload  distribution.  The  method  of  computing  air¬ 
loads  developed  in  this  investigation  makes  use  of  high-sneed  digital  computa¬ 
tion  to  retain  much  of  the  detail  of  the  vraks, 

rieasured  motion  and  strain  data  were  i  d  in  an  effort  to  separate  the 
airloads  problem  from  the  resnonse  problem  and,  thereby,  permit  concentration 
of  effort  in  this  investigation  on  the  formulation  of  an  adequate  aerodynamic 
representation.  It  was  found  ihat  separating  the  airloads  and  response  in 
this  manner  was  relatively  unsatisfactoiy  because  known  constraints  are  not 
satisfied  (e.g,,  the  first  harmonic  moment  about  the  flap  hinge  should  be  zero). 
If  the  motions  had  been  computed  with  the  airloads  instead  of  using  the  measured 
motions  and  separating  the  problem,  these  constraints  would  have  been  satis¬ 
fied. 


5 


The  method  developed  was  used  to  compute  the  airload  distribution  for 
three  flight  conditions  of  an  NASA  wind  tunnel  model  rotor  (described  in 
Refs,  1,  2  and  3)*  four  flight  conditions  of  the  Bell  HB-IA  full-scale  heli¬ 
copter,  and  two  flight  conditions  of  the  Sikorsky  H-31;  full-scale  helicopter. 
The  comparigpns  of  some  of  these  ccinputed  airloads  with  the  corresponding 
measured  airloads  are  presented  and  discussed  in  this  report. 

Sensitivity  of  the  computed  airloads  to  several  of  the  parameters  of 
the  computation  were  investigated  computationally,  and  the  results  are  pre¬ 


sented  and  discussed 


DEvSCRIPTION  OF  THE  METHOD 


The  Computational  Model 

It  is  assumed  that  the  rotor  is  operating  in  steady-state  flight. 

The  generated  wake  and  the  airloads  are,  therefore,  the  same  for  each 
revolution  of  the  blades,  i.e.,  the  form  of  the  vrake  is  periodic  and  the 
blades  ivill  "see"  the  same  wake  each  time  they  are  in  the  same  azimuth 
position.  It  is  also  assumed  that  the  effects  of  viscous  dissipation  in 
the  wake  can  be  neglected. 

Each  blade  of  the  rotor  is  represented  by  a  segmented  lifting-line 
(bound  vortex)  located  along  the  steacfy  deflected  position  of  the  quarter- 
chord.  The  number  of  segments,  fv  ,  and  the  length  of  each  segment  are 
arbitraryj  they  are  each  straight  and  of  constant  vortex  strength.  The 
lifting  line  is  considered  to  advance  in  a  stepwise  manner  through  ”  N  " 
equally-spaced  azimuth  positions. 

In  the  wake,  the  shed  and  trailing  vorticity  distributions  of  each 
blade  are  represented  by  a  mesh  of  segmented  vortex  filaments;  each  segment 
is  straight  and  of  constant  vortex  strength.  The  segmented  trailing  vortex 
filaments  emanate  from  each  of  the  end  points  of  the  lifting-line  segments. 
The  segmented  shed  vortex  filaments  intersect  the  trailing  filaments  in  a 
manner  such  that  the  end  points  of  both  are  coincident  (Fig.  l). 
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The  strengths,  d'/iv  ,  of  the  shed  elements  are  equal  to  the  change 
in  strength  of  the  bound  vortex  segments  between  successive  azimuth  stations 
and  are  deposited  in  the  flow  at  each  azimuth  station  of  the  bound  vortex. 

The  strengths,  d  /d'o  ,  of  the  trailing  vortex  elements  are  equal  to  the 
differences  in  strengths  of  adjacent  bound  vortex  segments  and  are  deposited 
in  the  flow  in  a  manner  such  that  they  connect  the  bound  vortex  end  points 
to  the  shed  vortex  end  points.  The  displacement  time  history  of  the  wake 
elements  is  specified,  but  this  specification  is  relatively  unrestricted. 

Thus,  in  this  method,  any  physically  realistic  distortion  of  the  wake  can 
be  incorporated  into  the  computation  of  the  airloads. 

The  Simultaneous  Equations 

For  each  azimuth  position,  the  airloads  are  computed  at  the  mid-points 
of  the  lifting-line  segments.  Thus,  the  airloads  are  computed  at  (  A  )•(  M  ) 
points  of  the  rotor  disk.  These  points  ^  of  the  disk  are  labeled  as 
shown  in  the  example  of  Fig.  2j  subscript ,  refers  to  position  in  the  disk. 
The  expression  for  the  lift  pe^unit  span,  for  a  blade  section  at  ar^’’  point,  ^ 
of  the  disk,  is 

K  -  " 

where 
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The  quantity  in  parentheses  is  the  instantaneous  "effective”  angle-of -attack 
at  the  three-quarter-chord  point,  where  is  the  instantaneous  geometric 
angle-of -attack  with  respect  to  the  tip-path-plane,  the  "effective" 

plunging  velocity  with  respect  to  the  tip-path-plane,  and  UTj^  the  velocity 
coii?)onent  normal  to  the  tip-path-plane  induced  by  all  the  wake  vorticity  and 
the  bound  vorticity  of  other  blades.  The  strengths  of  the  bound  vortex 
elements  are  considered  to  be  the  unknowns  for  the  method  of  solution  developed. 

The  velocity  coirponent  in  a  given  direction  induced  at  a  point  by  an 
arbitrarily  oriented  straight  vortex  filament  of  constant  strength  is  given 
by  the  Biot-Savart  law  as: 

(3) 

where  Y  is  the  constant  vortex  strength  of  the  filament,  and  the  coefficient 
^  is  a  function  only  of  the  coordinates  of  the  point  where  the  velocity  is 
being  corrpited  and  the  coordinates  of  the  voilex  filament  end  points.  The 
velocity.  UJT  .of  Eq.  2  can  be  coiirouted  by  summing  the  contributions  (given 
by  Eq.  3)  of  each  individual  vortex  element  of  the  blades  and  wake  as  indicated 
by  Eq.  U. 

(w 

However,  each  wake  vortex  element  strength  Tf  of  Eq.  k  is  just  the  difference 
in  strength  of  two  of  the  unknown  bound  vortex  segment  strengths,  (i.e., 

Y  ■  ,  idiere  <V-  and  h  are  two  values  of  the  disk  position 

subscript,  ).  These  Y  relationships  are  depicted  in  Fig.  1  relative  to  a 
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point  1|^  of  the  rotor  disk.  If  these  differences  are  substituted  for  the 
Y  ,  Eq,  k  can  be  expanded  and  the  unknown  grouped  and  factored  out* 
Thus,  the  can  be  expressed  as 


nN 


In  Eq.  the  are  the  rinknown  bound  vortex  strengths  and  each  is  the 

sum  of  all  the  coefficients  of  a  common  T;^  in  Eq.  U.  substituting 
the  given  by  Eq.  ^  in  Eq.  2  and  letting 

following  expression  for  is  obtained: 

A  set  of  linear  simultaneous  nonhoraogeneous  equations  with  constant  coefficients 
is  obtained  when  Eq.  6  is  written  for  each  point,  ,  of  the  rotor  disk  at 
which  the  airload  is  to  be  con^juted.  Thus,  the  airloads  problem  is  reduced 
to  obtaining  the  solution  to  this  set  of  equations  for  which  the  unknowns  are 
the  (  v\  )»■(  N  )  bound  vortex  strengths  • 

Solution  of  the  Equations 

An  iterative  procedure  idiich  is  described  below  is  used  to  solve  the 
set  of  simultaneous  equations.  First,  Eq,  6,  the  general  equation  of  the  set, 
is  rearranged  and  put  in  the  following  form  with  the  unknowns  on  the  right  side  , 
of  the  equation: 
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v\H  ^ 

where 

fu,  *  l.O 

Thus,  the  set  of  equations  to  be  solved  is: 


1,--  (1.011:  .f 

^^Ti.  +  (r_T,Tj  i- . 

!,  •-  + 

*■  1  * 

4  .  .  .  . 

*  •  •  *  • 

K 

‘  T 

I.-  i- 

4-  (,l‘0)"\\4  •  .  *  . 

•  •  +^x,»,.nT,.h 

in.N'  + 


•  •  -  •  +O.0)T;1.n 


For  the  first  iteration  only,  the  set  of  equations  is  assumed  to  be  con?)letely 
diagonalized  (i.e.,  the  off-diagonal  (f  »  are  assumed  equal  to  zero);  thus, 
from  the  first  iteration,  the  are  equal  to  the  ■vdiere  the  siperscript 
represents  the  iteration  number.  The  second  and  succeeding  approximations  to 
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the  Tl  are  obtained  by  an  iterative  procedure  governed  by  the  following 


equation  which  is  solved  in  succession  for  with  I 


The  following  is  a  brief  word  description  of  the  iterative  procedure  governed 
by  Eq.  9: 


Beginning  with  the  first  equation  in  the  set,  each  is  solved  in  succession 
for  the  unknown  T  on  the  diagonal  by  using  for  all  T  »  to  the  left  of  the 
diagonal  the  approximations  to  )  thus  far  obtained  in  the  present  (j  ) 

/  I  V 

iteration  (i.e.,  the  where  k  represents  the  equation  number) |  and 

for  all  the  I  \  to  the  right  of  the  diagonal,  the  approximate  values  obtained 
in  previous  )  iteration  (i.e.,  the  ).  In  Eq.  9,  the  terms  to 

the  left  of  the  diagonal  are  represented  by  the  first  summation  and  those  to 
the  right  by  the  second  summation. 


Discussion 

The  following  are  the  major  assumptions  which  have  been  made  in  developing 
the  computational  model  and  procedure  and  have  made  it  possible  to  reduce  the 
airloads  problem  to  a  set  of  simultaneous  equations; 

1,  The  rotor  has  been  in  steady  flight  long  enough  for  the  wake  to 
become  periodic,  i.e.,  the  generated  wake  is  the  same  for  each 
revolution  of  the  blades. 
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2,  The  effects  of  viscous  dissipation  on  wake  vortex  strengths  can 
be  neglected. 

3.  The  blades  can  be  represented  by  lifting  lines. 

U.  The  trailing  and  shed  vorticity  distributions  of  the  wake  can  be 
represented  by  a  mesh  of  segmented  vortex  filaments. 

5.  The  displacement  time  history  of  these  wake  elements  is  specified. 

The  framework  of  the  method  has  been  made  sufficiently  general  to  permit 
incorporation  of  the  effects  of  the  following: 

1.  Number  of  blados 

2.  Finite  span 

3.  Planfom  taper 
U.  Twist 

Variable  lift  curve  slope 

a,  Mach  number 

b.  Reynolds  number 

6.  Stall  and  reversed  flow 

7.  Blade  elastic  deformations 

8.  Control  system  deformation 

9.  Blade  rigid  boc^y  motions 

10.  Mutual  interference  of  rotors 
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stall  effects  are  introduced  into  the  computations  during  the  solution 
of  the  simultaneous  equations  by  determining  the  effective  angle  of  attack 
and  limiting  the  accordingly.  For  the  computations  which  were  made,  the 

lift  coefficient  was  assumed  to  be  constant  and  equal  to  its  value  at  the 
stall  angle  for  angles  of  attack  greater  than  the  stall  angle.  The  limiting 
"^1^  (i.e.,  the  stalled  value  of  )  for  these  calculations  were  then 

computed  from  Eq.  6  by  substituting  the  airfoil  section  stalling  angle  of 
attack  for  the  effective  angle  of  attack,  j  thus, 

•P 

the  limiting  4  were  computed  from  the  following: 

t-- 

It  is,  however,  possible  within  the  present  framework  of  the  method  to  use 
more  representative  lift  coefficient  variations  with  angle  of  attack  than 
the  one  chosen  here.  It  should  be  noted  that  incorporating  stall  in  this 
manner  assumes  that  all  of  the  stalled  section  lift  is  circulatory  and  there 
is,  therefore,  shed  and  trailing  vortex  elements  in  the  wake  whose  strengths 
are  determined  by  these  stalled  values  of  • 

In  reality,  when  an  airfoil  section  is  stalled,  the  predominant  part  of 
the  airload  on  the  section  may  be  noncirculatory  and  due  to  the  separated  flow. 
During  the  stalled  condition,  the  vorticity  introduced  in  the  wake  may  be 
relatively  small  and  of  a  considerably  different  nature  than  that  of  an  vui- 
stalled  section.  Some  indication  of  the  magnitude  of  the  separated  flow  forces 


experienced  is  evident  when  the  total  measured  airloads  are  examined.  If 
the  local  dynamic  pressure  and  measured  airload  are  used  to  compute  the 
corresponding  normal  force  coefficient,  it  is  found  to  be  extremely  large 
for  some  points  of  the  rotor  disk.  These  large  "measured"  normal  force 
coefficients  occur  in  regions  of  the  rotor  disk  where  the  blades  are  supposedly 
stalled.  These  forces  are  analogous  to  the  normal  force  on  a  flat  plate  at 
very  large  angles  of  attack. 

Reversed  flow  effects  are  accounted  for  in  the  cor^jutation  approximately 
in  that  the  sign  of  the  lift  force  and  the  sense  of  the  bound,  shed  and  trail¬ 
ing  vorticity  are  always  consistent  with  the  local  tangential  velocity,  V  , 
and  the  "effective"  relative  angle-of -attack  of  the  section  as  it  approaches 
and  passes  through  the  region  of  reversed  flor^.  However,  for  the  computations 
which  were  performed,  the  wake  remained  attached  to  the  bound  vortex  of  the 
blade  (although  its  position  and  orientation  corresponded  to  the  local  inflovr 
velocity)  and  the  bound  vortex  remained  at  the  quarter-chord  while  the  blade 
approached  and  passed  through  the  region  of  reversed  flow.  In  reality,  when 
a  section  is  operating  unstalled  in  reversed  flow,  the  center  of  pressure  is 
approximately  at  the  three-quarter-chord  and  the  bound  vortex  should  be  placed 
at  that  station  and  the  Induced  velocity  calculated  at  the  quarter-chord, 

Mach  number  (less  than  one)  and  Reynolds  number  effects  can  be  accounted 
for  as  variations  of  the  lift  curve  slope  of  the  local  blade  sections.  For  a 
given  rotor  and  operating  conditions,  the  Ifech  number  and  Reynolds  number  distri¬ 
butions  over  the  rotor  disk  are  fixed  (known) j  thus,  these  effects  can  be 
incorporated  into  the  computations  by  introducing  the  proper  distribution  of 
lift  curve  slope  over  the  disk. 


Parameters  of  the  Computations 


A  computational  investigation  was  made  to  determine  the  sensitivitj'"  of 
the  computed  airloads  to  several  "parameters"  of  the  method;  They  are, 
the  number  of  revolutions  of  wake  used,  the  number  of  radial  segments  used 
to  represent  the  blades,  and  the  distortion  of  the  wake.  This  computational 
investigation  was  made  for  the  NASA  model  rotor  at  an  advance  ratio  of  0,l5, 
Ejcperimental  azimuthal  load  distributions  are  shown  on  Figure  3  for  this  case. 

It  was  found  for  this  rotor  and  operating  condition  that  truncating 
the  wake  at  three  revolutions  put  the  wake  termination  approximately  one- 
and-one-half  rotor  diameters  away  from  the  rotor  disk  and  accounted  for 
practically  all  of  the  wake-induced  velocities j  extending  the  wake  to  five 
revolutions  only  resulted  in  a  maximum  change  of  two  percent  in  the  magnitudes 
of  the  lof/rest  harmonics  of  the  computed  airloads.  This  v/as  expected  because 
the  additional  wake  is  relatively  far  from  the  rotor  and  the  induced  velocities 
due  to  it  are  relatively  uniform  over  the  diskj  whereas,  the  induced  velocity 
distribution  due  to  the  nearer  or  most  recent  wake  has  a  considerable  variation 
over  the  rotor  disk  and  is  the  nrimarj^  source  of  the  higher  harmonic  components 
of  the  airloads.  It  is  believed  that  the  amount  of  wake  (expressed  in  terms 
of  a  corresponding  number  of  rotor  revolutions)  which  is  necessary  to  adequately 
predict  the  induced  velocities  at  the  blades  will  depend  on  the  given  rotor  and 
its  flight  condition. 

For  this  same  rotor  and  operating  condition,  the  sensitivity  of  the 
computed  airloads  to  the  number  of  radial  segments  used  to  represent  the  blades 
vras  investigated  by  computing  the  airloads  for  both  five  and  nine  radial 
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segments.  A  comparison  of  the  results  of  these  computations  is  presented  in 
Figure  Ij..  At  this  time  it  is  not  nossible  to  draw  any  specific  conclusions 
from  these  results  other  than  to  note  that  a  change  in  the  number  of  radial 
segments  from  five  to  nine  for  this  rotor  and  operating  condition  can  have 
some  effect  on  the  airloads. 

For  both  the  five  and  nine  radial  segment  computations,  the  same  first 
harmonic  .flanning  motions  were  used.  However,  the  constraint  of  zero -moment 
about  the  flapping  hinge  due  to  the  first  harmonic  component  of  the  airload 
is  not  satisfied  for  either  of  these  cases;  this  constraint  arises  from  the 
flapping  equation  of  motion.  For  a  proper  comparison  of  the  airloads  computed 
Tfri-th  five  and  nine  radial  segments,  this  flapping  moment  constraint  should  be 
applied  in  each  case.  This  could  be  done  by  including  the  flapping  equations 
of  motion  in  the  set  of  equations  and  solving  simultaneously  for  the  airloads 
and  the  flapping  motion. 

Sensitivitir  of  the  computed  airloads  to  distortion  of  the  wake  was  also 
investigated  for  this  rotor  and  operating  condition  by  using  several  relatively 
arbitrary  distortions.  The  simplest  wake  distortion  used  corresponded  to  a 
radial  distribution  of  the  zeroth  harmonic  cotnponent  of  the  inflow 
velocity  distribution  in  the  nonrotating  system.  The  most  elaborate  corres¬ 
ponded  to  the  use  of  harmonics  of  this  velocity  distribution  through  the  fourth 
with  magnitudes  up  to  thirty  percent  of  the  zeroth  harmonic.  Common  to  each 
of  these  wake  distortions  was  the  assumption  that  the  velocity  of  transport 
of  each  vrake  element  end  point  was  constant  with  time.  The  velocity  used 
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for  each  segment  end  point  was  taken  to  be  the  vector  sura  of  the  rotor  trans¬ 
lational  velocity  and  the  component  of  inflow  velocity  normal  to  the  rotor 
disk  (tip-path-plane)  at  the  point  of  the  disk  where  the  segment  end  point 


originated  (i.e.,  the  end  points  were  assumed  to  retain  their  initial  velo¬ 
cities).  Thus,  at  any  instant  of  time,  the  vjake  element  end  points  defined 
a  distorted  helical  surface.  The  wake  of  each  blade,  at  any  instant  of  time, 
is  different  and  changes  with  azimuth  position  to  correspond  to  the  instanta¬ 
neous  wake  configuration  for  that  azimuth  position. 

It  was  found  that  the  computed  airloads  were  sensitive  to  these  distor¬ 
tions  of  the  wake  for  this  flight  condition  of  this  rotor.  It  should  be  noted 
here,  again,  that  the  distortions  of  the  vialce  (i,e,,  displacement  time  history 
of  the  wake)  which  can  be  used  in  this  method  of  computing  rotor  airloads 
are  relatively  unrestricted.  There  remains,  however,  the  question  as  to  what 
displacement  time  history  should  be  used  for  a  particular  rotor  and  operating 
condition.  Using  a  uniform  inflow  distribution,  determined  from  momentum  con¬ 
siderations,  to  define  the  time  history,  has  resulted  in  computed  airloads 
which  agree  quite  vrell  with  the  measured  airloads.  This  is  evidently  a  good 
first  approximation  for  the  conditions  analyzed j  therefore,  the  addition  of 
a  few  of  the  lower  harmonic  variations  of  the  inflow  distribution  (in  the 
nonrotating  system)  may  be  all  that  is  necessary  to  adequately  define  the 
wake  displacement  time  history.  It  is  possible,  by  making  a  fex'i  modifications 
to  the  present  method,  to  compute  in  addition  the  periodic  induced  velocity 
distribution  in  the  nonrotating  system. 
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The  method  as  developed  for  computing  the  airload  distribution  has  been 
programmed  (in  FORTRAl'I  for  the  IBH  7Ch)  such  that  it  is  relatively  simple  to 
use.  The  printed  output  obtained  from  the  program  includes  the  following 
distributions  given  as  functions  of  the  radial  and  azimuthal  position  in  the 
disk, 

1,  Airload  on  the  blades 

2,  Strength  of  the  bound  vorticity 

3,  Effective  angle  of  attack 
li.  Stalled  angles  of  attack 

Induced  velocitj'-  at  blade 
6,  Induced  drag 

Items  1,  2,  5  and  6  are  also  harmonically  analyzed  and  the  Fourier  sine  and 
cosine  coefficients  printed  out. 
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mSURED  RESULTS  USED  FOR  CORRELATION 


Source 

The  measured  airloads  and  motions  which  were  used  to  test  the  adequacy 
of  the  computational  model  and  procedures  were  selected  from  three  sources: 
full-scale  flight  tests  of  the  HU-IA  helicopter  by  Bell  Helicopter  Corporation 
(Ref.  36),  vdnd  tunnel  tests  of  a  model  rotor  by  the  NASA  (P^fs.  1  -  3)j  and 
full-scale  flight  tests  of  the  H-3U  helicopter  by  the  NASA  (Ref,  37 )•  The 
HU-IA  has  a  two-blade  teetering  rotor  and  the  H-3l:  a  four-blade  fully  articu¬ 
lated  rotor j  both  have  constant-chord  twisted  blades.  The  NASA  model,  which 
is  a  two-blade  teetering  system,  has  constant- chord,  nontwisted  blades.  The 
disk  loadings  for  these  three  rotors  are  approximately  as  follows:  HU-IA, 

U.O  psf j  H-3U,  U,9  psf j  NASA  model,  -2.3  psf. 

Blade  Motion  Analysis 

The  following  measured  quantities  relating  to  the  blade  motions  were 
used  in  this  analysis; 

1.  Blade  root  rotation  about  flapping  hinge  as  function  of  azimuth, 

2.  Blade  root  rotation  about  pitch  axis  as  function  of  azimuth, 

3.  Flapwise  bending  moments  along  blade  span  as  function  of  azimuth, 

U.  Torsion  moments  along  span  as  function  of  azimuth. 

Items  1  and  3  were  harmonically  analyzed. 
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The  blade  motions,  as  used  in  the  present  approach  for  computing  the 
airloads,  are  classified  as  either  plunging  or  rotational.  The  plunging 
displacements  are  defined  to  be  the  normal  displacements  of  the  quarter- 
chord  with  respect  to  the  tip-path-plane.  The  rotational  displacements 
are  about  the  quarter- chord  and  with  respect  to  the  tip-path -plane. 

The  radial  distributions  of  the  plunging  displacements  at  each  harmonic 
of  the  rotational  speed,  sine  and  cosine,  are  computed  by  double  integration 
of  a  Fourier  series  fit  to  the  radial  curvature  distributions;  the  curvature 
distributions  are  obtained  from  the  bending  stiffness  and  measured  moment 
distributions  at  each  haimonic,( obtained  from  harmonic  analysis  of  measured 
moments).  The  constants  of  integration  are  evaluated  from  the  end  conditions 
which  include  the  measured  flapwise  root  rotation  at  each  harmonic.  The  slopes 
and  plunging  velocities  of  the  blades,  x^hich  are  used  directly  in  the  airloads 
computation,  are  obtained  by  differentiating  the  displacement  distribution 
Tilth  respect  to  the  radius  and  time,  respectively.  The  first  harmonic  flapping, 
sine  and  cosine,  is  obtained  ty  performing  an  orthogonal  analysis  of  the  radial 
distribution  of  the  first  harmonic  plunging  displacement;  the  orthogonal  fiuictions 
used  are  the  first  five  bending  mode  shapes  of  the  blade.  This  procedure  for 
computing  the  blade  deflections  from  the  measured  moment  distributions  is,  in 
principle,  similar  to  that  reported  by  Sikorsky  (Ref,  38). 
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For  both  the  rlU-lA  and  the  the  only  T-easured  data  available  to  • 

determine  the  pitching  rotational  motions  of  the  blade  sections  are  the 
pitch  root  rotation  and  the  torsional  moments  at  only  two  radial  stations 
(15  and  50  percent  of  radius) .  At  each  azimuth  station,  tvro  simultaneous 
equations  are  set  up  and  solved  for  the  tip  deflection  in  each  of  the  first 
two  normalized  (for  unit  tip  deflection)  cantilever  torsional  modes  by  malcing 
use  of  their  moment  distributions  and  the  measured  torsional  moments  at  each 
of  the  two  radii.  The  rotational  deflections  are  then  computed  from  these 
tip  deflections  in  each  torsion  mode,  the  computed  torsional  mode  shapes, 
and  the  nitch  root  rotations. 

In  using  the  above  procedures  for  comnuting  the  blade  motions  from  the 
measured  root  displacements  and  the  bending  and  torsion  moment  distributions, 
the  effects  of  intermodal  counling  on  the  motions  have  been  neglectedj  the 
effects  ma3'',  however,  be  imnortant.  If  the  effects  of  intermodal  coupling 
were  to  be  included,  it  would  be  necessary  to  increase  the  accuracy  of  the 
motion  data  considerably  and  perform  more  extensive  calibrations. 

It  should  be  noted  here,  again,  that  this  approach  of  using  the  measured 
motion  and  strain  data  to  compute  the  airloads  was  taken  onljr  as  an  attempt 
to  separate  the  airloads  problem  from  the  over-all  aeroelastic  problem  and, 
thereb;/,  enable  the  focusing  of  effort  on  the  formulation  of  an  adequate  aero¬ 
dynamic  representation.  Although  some  knovm  constraints  were  not  satisfied 
because  of  this  apnroach,  it  was  successful  in  that  the  objective  was  attained 
now  that  the  basic  aerodynamic  representation  has  been  developed,  effort  can 
be  spent  on  refinements  and  the  complete  aeroelastic  problem. 
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RESULTS 


Measured  and  computed  airloads  are  presented  for  the  following 
configurations  and  flight  conditions: - ■ - - . 

MSA  Model]  fA  =  0.08,  0.15,  0.29 
w  -  0.18,  0.29 

HU-U]  j.K  =  0.21,  0.26 

Tiio  of  these  viill  be  eit^Dhasized  —  the  NASA  model  at  yt-*  =  0,1^  and  the  HU-IA 
at  ytv.  =  0.26.  For  these  two  cases,  the  airloads  are  compared  in  two  ways; 

(1)  azimuthal  distributions  at  each  radial  station,  and  (2)  radial  distributions 
of  the  steadj'-  and  the  sine  and  cosine  components  of  each  of  the  first  eleven 
haimonics.  The  two  forms  of  presentation  are  used  to  enrohasize  the  fact  that 
apparent  good  agreement  with  respect  to  azimuthal  distributions  does  not 
necessarily  guarantee  good  agreement  with  respect  to  particular  harmonics. 
Comparisons  should  always  be  made  with  respect  to  the  spamdse  distributions 
of  the  harmonics  of  the  airloads  because  it  is  believed  that  this  is  the  form 
most  suitable  for  the  derivation  of  the  forces  transmitted  to  the  fuselage  and 
the  generalized  forces  for  estimating  oscillatory  blade  stresses. 

All  graphs  of  oscillatorjr  airloads  for  a  given  rotor  and  flight  condition 
are  presented  on  the  same  scale  to  facilitate  inspection  of  relative  amplitudes 
at  the  various  harmonics  and  advance  ratios.  Curves  have  been  passed  through 
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the  experimental  and  theoretical  points  to  indicate  the  form  of  the  spam-rise 
distribution  of  airload  components.  Unfortunately,  the  points  are  not  always 
sufficiently  dense  to  define  a  unique  curve. 

Mach  number  and  Reynolds  number  corrections  have  not  been  made,  although 
these  could  be  introduced  without  modifying  the  calculation  procedure  developed. 

A  constant  lift  curve  slope  of  ?.73  per  radian  was  used,  except  in  those  calcu¬ 
lations  in  which  stall  effects  were  included.  Stall  effects  xrere  included  in 
the  calculations  made  for  the  KU-IA  at  yA  =  0.26  by  limiting  the  circulatory 
lift,  as  previously  described, 

A  typical  induced  velocity  distribution,  obtained  from  one  of  the  compu¬ 
tations  (NASA  model,  yU,  =  0.15,  nine  radial  segments),  is  presented  in  Table  1. 
The  rox^s  in  the  table  are  the  radial  distributions  from  the  root  to  the  tip  of 
the  induced  velocity  at  each  azimuth  position  from  4^  =  0°  to  =  360°,  in 
15°  increments.  These  are  the  induced  velocities  on  the  blade  as  it  rotates 
and  are,  therefore,  in  the  rotating  coordinate  sj'-stem. 

It  is  interesting  to  note  that  the  computed  induced  velocity  distributions 
are  alxrays  such  as  to  oppose  the  retreating  blade  tip  stall  which  has  always 
been  predicted  on  the  basis  of  an  assximed  uniform  inflow  distributionj  this  is 
more  clearly  evident  from  Table  2  in  x^hich  the  distribution  of  the  ratios  of 
the  instantaneous  effective  angle  of  attack  to  an  angle  of  ll;°  is  presented. 

Also  to  be  noted  are  the  positive  induced  velocities  (up-wash)  and  their  relative 
position  in  the  disk.  Regarding  this  point,  it  should  be  remembered  that  these 


are  the  induced  velocities  on  the  blade  for  its  various  positions  in  the 
rotor  diskj  furthermore,  the  up-wash  through  the  disk  xirhich  has  been  observed 
and  photographed  by  other  investigators  is  the  result  of  the  time-averaged 
induced  velocities  at  fixed  points  of  the  disk  as  the  blades  rotate  (i.e., 
the  steady  component  of  the  periodic  induced  velocities  in  the  nonrotating 
coordinate  -sy-stemX* 

NASA  Model,  =0.1$ 

The  wake  time  history  used  for  these  computations  was  defined  by  the 
tunnel  speed,  tip-path-plane  inclination,  and  an  assumed  uniform  value  of 
induced  velocity  based  on  momentxun  theory.  The  wake  used  corresponded  to 
three  revolutions  of  the  rotor.  The  measured  collective  pitch  angle  was  used. 
For  the  computation  using  five  radial  control  points,  the  first  harmonic 
flapping  was  estimated  in  order  to  satisfy,  approximately,  the  condition  that 
the  first  harmonic  airload  moment  about  the  teetering  pin  be  zero.  Only  the 
first  harmonic  flapping  motions  were  included  — -  the  blade  elastic  deformations, 
pitching  and  higher  harmonic  flapping  motions  were  neglected. 

The  computed  and  measured  azimuthal  airload  distributions  are  presented 
in  Fig.  3»  The  computed  curves  are  based  on  five  radial  control  points  and, 
to  proxTide  a  more  direct  comparison  of  the  oscillatory  corponents,  they  are 
corrected  for  the  difference  between  the  computed  and  experimental  steacty’ 
airloads.  Good  over-all  agreement  is  indicated,  A  much  more  stringent 


comparison  is  made  in  Fig.  U  vihere  the  radial  distributions  of  the  sine  and 
cosine  components  at  each  harmonic  of  the  airloads  are  presented. 

Calculated  and  measured  steady  airload  distributions  are  shown  on  Fig. 
l;a.  The  computed  curve  based  on  five  radial  control  points  (in  fact,  the 
same  radial  stations  at  v/hich  the  airloads  were  measured)  is  showm  as  a 
dashed  line,  and  the  curve  based  on  nine  radial  points  is  shoivn  as  a  dotted 
linej  a  solid  line  is  draTm  through  the  experimental  data  in  these  graphs. 

The  form  of  the  steady  airload  distribution  is  obtained,  but  the  computed 
magnitudes  are  too  large. 

First  harmonic  cosine  and  sine  load  distributions  are  shorn  on  Fig.  l;b. 
Attention  is  again  called  to  the  fact  that  the  first  harmonic  flapping  was 
chosen  to  yield  a  cosine  airload  distribution  having  a  nearly  zero  moment  about 
the  flapping  hinge  for  the  computation  based  on  five  radial  control  points. 

It  is  immediately  evident  that  the  computed  curve  based  on  nine  radial  control 
points  and  this  same  assvuned  flapping  does  not  satisfy  the  cosine  lift  moment 
constraintj  however,  the  experimental  curve  drawn  does  not  satisfy  this  con¬ 
straint  either.  In  Ref.  2,  the  expeid-mental  load  distribution  was  evidently 
made  to  satisfy  the  lift  moment  constraint’  by  requiring  sizable  negative  lifts 
at  r/j^  “  both  sine  and  cosine  components.  Neither  of  the  computed  cosine 

components  agrees  with  the  measured  distribution  as  drawn.  The  conpited  first 
harmonic  sine  components  of  airload  obviously  do  not  satisfy  the  lift  moment 
constraint,  and  improvement  in  agreement  would  be  obtained  by  ir5>osing  this 
condition.  The  computed  first  harmonic  airloads  were  found  to  be  quite  sensitive 
to  the  first  harmonic  flapping  angle. 
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Second  harnonic  components  of  the  airloads  are  sh'own  on  Figure  l;c« 

Theory  and  experiment  agree  as  to  distribution  and  reasonably  well  with 
respect  to  magnitude. 

Third  harmonic  components  are  shown  on  Figure  lid.  Reasonable  agreement 
was  obtained.  There  is  an  interesting  point  here.  Although  the  first  un¬ 
coupled  uns;/mraetric  bending  mode  frequency  of  the  model  is  nearly  three  per 
rev  (Ref.  2),  the  calculated  airload  distribution,  in  which  the  elastic 
motions  vjere  ignored,  agrees  reasonably  x-rell  with  the  measured  values,  VJhether 
or  not  this  is  the  result  of  compensating  errors  is  not  known  at  this  time. 

Figures  Ije  through  are  graphs  of  the  congsuted  and  measured  airloads 
of  the  fourth  through  the  eleventh  harmonic.  In  nearly  all  cases  the  measured 
and  calculated  values  agree  reasonably  well.  Because  only  first  harmonic  blade 
motions  were  included,  the  computed  airloads  at  harmonics  beyond  the  third 
are  entirely  due  to  the  induced  velocity  distribdtion.  The  experimental  air¬ 
loads  for  harmonics  beyond  the  sixth  v;ere  obtained  by  harmonic  analysis  of 
the  azimuthal  distributions  given  in  Reference  Ij  no  correction  was  made  for 
instriimentation  response,  A  previous  report  related  to  this  rotor  (Ref,  3) 
indicates  that  the  instrumentation  correction  is  trivial  at  frequencies 
corresponding  to  the  sixth  harmonic,  but  at  the  eleventh  harmonic  the  measured 
amplitudes  would  be  approximately  20/»  too  hi^. 

In  summary,  the  calculated  and  measured  airload  distributions  for  the 
model  rotor  showed  qualitative  agreement,  A  relatively  simple  wake  time 
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histoiy  was  usedj  blade  flexibility  was  neglected,  and  only  the  first  harmonic 
rigid-bodjr  flapping  motions  were  used.  Furthermore,  knovm  constraints  implicit 
in  satisfying  the  flapping  equation  of  motion  were  not  rigorously  satisfied. 


HU-IA  at  - - 

The  computed  and  measured  airloads  for  the  HU-IA  at  /-(.  =  0,26  are  pre¬ 
sented  in  Figures  5  and  6.  The  wake  time  history  used  for  these  computations 
was  defined  by  the  free  stream  velocity,  tip-path-plane  inclination,  and  an 
assumed  uniform  induced  velocity  derived  from  momentum  considerations,  A 
wake  corresponding  to  two  rotor  revolutions  was  used.  Measured  values  of 
blade  root  rotation  and  blade  strains  were  used  to  determine  blade  positions, 
shapes,  and  velocities;  however,  the  effects  of  intermodal  coupling  were 
neglected  when  determining  these  motions.  The  results  of  two  sets  of  calcu¬ 
lations  are  presented  in  Fig,  6:  the  first,  in  which  only  the  steady  deflec¬ 
tions  and  rigid-body  first  harmonic  motions  were  included;  and  the  second, 
in  which  all  the  rigid  and  elastic  motions  were  included.  For  the  case  in 
which  the  higher  harmonic  motions  were  inclu.ded,  the  motion  and  strain  data 
traces  were  inspected  and  those  motions  and  strains  VThich  were  represented  by 
trace  deflections  of  less  than  0,02  inches  were  discarded.  Stall  effects  vrere 
included  in  the  manner  discussed  previously.  Airloads  were  conqjuted  at  eight 
radial  stations,  five  of  which  are  the  stations  at  which  the  airloads  were 
measured. 


28 


The  azimuthal  variations  of  the  measured  and  computed  airloads  at  each 
radial  station  are  presented  in  Fip.  5?  the  computed  airloads  presented  here 
are  from- the  computation  which  included  all  the  rigid  and  elastic  motions. 
Fieasonable  agreement  is  obtained. 

The  steady  component  of  the  airloads  are  sho-'.'m  on  Fig,  6a,  and  it  can 
be  seen  that  the  measured  and  calculated  values  are  in  good  agreement. 

First  harmonic  airloads  are  shoim  on  Fig.  6b.  Although  the  distri¬ 
bution  of  the  computed  first  harmonic  sine  component  agrees  quite  well  id-th 
the  measured  distribution,  the  distri.bution  of  the  computed  cosine  component 
does  not,  especially  over  the  inboard  position  of  the  blade.  Including  the 
elastic  motions  and  the  harmonics  of  flapping  above  the  first  had  an  effect 
on  the  computed  airloads,  but  there  are  no  clear  grounds  for  assessing  the 
necessity  for  including  elastic  deformations. 

Figure  6c  is  a  graph  of  the  second  harmonic  airloads.  The  computed 
cosine  component  agrees  reasonably  vjell  vri.th  the  measured  cosine  component 
but  the  sine  component  agreement  is  poor.  Again,  elastic  deformations  used 
had  relatively  little  effect. 

Third  harmonic  airloads  are  shovm  on  Fig.  6d.  Here,  there  is  a  temptation 
to  argue  that  inclusion  of  the  elastic  motions  and  higher  harmonics  of  flapping 
is  significant  —  especially  because  there  is  a  natural  flapid.se  bending 
frequency  in  the  vicinity  of  the  third  harmonic  and,  therefore,  the  third 
haimonic  bending  motions  are  relatively  large.  However,  this  argument  cannot 
be  sustained  on  the  basis  of  the  airloads  at  a  single  harmonic  because  the 
motions  at  a  given  harmonic  not  only  contribule  directly  to  the  airloads  at 
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that  harmonic  but,  also,  contribute  c^irectly  to  the  airloads  at  the  harmonic 


above  and  belo^r  the  given  harmonic.  For  example,  examine  the  comparison  of 
the  calculated  and  measured  airloads  at  the  fourth  hamonic.  Fig,  6e  (the 
third  harmonic  blade  motions  contribute  directlj?'  to  the  fourth  harmonic). 

Here,  the  neglect  of  the  elastic  and  higher  harmonic  iTanning  notions  apparently 
jsLelds  better  agreement. 

The  fifth  and  higher  harmonics  are  shown  on  Figs.  6f  through  6j<  .  These 
are  nresented  as  a  natter  of  interest  in  that  they  are  relatively  small  in 
aimlitude.  The  significance  of  the  detailed  comparison  of  the  calculated  and 
measured  results  is  obscured  by  the  question  of  measurement  accuracy/, 

I-IU-ll  at  ^  =0,21 

Measured  and  calculated  harmonic  components  of  the  airloads  on  the 
•IU-1.A.  at  ^  =  0.21  are  shoi-m  in  Figure  7.  Stall  effects  were  judged  to  be 

small  and  were  neglected,  and  only  first  harmonic  rigid-body  motions  were 
retained  in  the  computations.  The  agreement  between  calculated  and  measured 
results  is  about  the  same  as  that  observed  for  the  previous  case  (ifu-lA  at 
M  =  0.26). 


NASA  Model;  M  =  0.08  and  0,29 

The  coiirouted  airloads  are  based  on  measured  collective  pitch,  shaft, 
and  first  harmonic  flapping  angles,  A  viake  corresponding  to  five  rotor 
revolutions  was  retained  for  the  =  0,08  case  and  two  revolutions  for  the 
M  =  0,29  case.  For  each  case  the  wake  was  positioned  on  the  assumption  of 
a  uniform  inflow  velocity  derived  from  momentum  considerations. 
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A  comparison  of  the  theoretical  and  exmeriraental  airloads  is  presented 
in  figures  3  and  9.  In  both  cases  the  mean  lood  was  overpredicted. 

It  is  evident  that  the  theoretical  first  harmonic  airload  distributions 
do  not  satisfy  the  requirement  that  the  moment  of  the  airload  about  the 
:nanning  hinge  be  zeroj  but  the  measured  cosine  airload  at  =  0.29  does 

not  satisf^r  this  requirement  either.  Introduction  of  the  rigid-body  f].anning 
equation  into  the  airload  analysis  would  force  the  theoretical  distribution 
to  satisfy  the  airload,  noaient  condition  but,  then,  the  calculated  f].aDping 
angles  would  not  agree  with  the  measured  values,  Resolution  of  this  apparent 
incompatibility  could  be  obtained  if  the  root  rotation  (which  is  an  item  that 
was  measured)  contained  an  apnreciable  -first  harnonic  bending  slope  suner- 
imnosed  on  the  rigid-body  flapning.  dloviever,  the  displacement  data  presented 
in  Ref.  2  indicates  this  is  unlikely. 

Relatively  good  agreement  was  obtained  for  tlie  second  harmonic  cosine 
components,  but  the  s-'-all  sine  coranonents  are  not  well  nredicted.  It  should 
be  noted  here,  again,  that  the  first  harmonic  flapning  velocity  and  c;;rclic 
pitch  influence  the  airloads  quite  strongly. 

There  is  reasonable  agreement  between  calculated  and  measured  values 
for  the  thdrd  harmonic.  For  this  particular  rotor,  the  third  hai*monic  is 
appreciably  larger  at  yA,  =  0,08  than  at  0,29. 

Agreement  between  measured  and  calculated  curves  is  good  for  the  .fourth 
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=  O.Ofj  ir;  undemredicbed 


and  fifth  harmonic  comnoncnts,  but  the  sixth  at 
at  the  outboard,  stations. 


Comparison  of  the  airloads  for  the  seventh  through  eleventh  lianionics 
is  also  nresented.  It  was  |■'ointed  out  ■previous Ip  that  those  ineasu.red  com¬ 
ponents  x-rere  obtained  bv  a  har -ionic  analpv^^'ie  of  the  azir.iuth.al  di.otrj'-bution 
of  the  airloads  fiven  granhically  in  hef.  1,  Because  n^'  corrections  were  -.■oade 
for  instrumentation  resnonse  (and  hef.  3  indicates  that  the  recorcied  eleventh 
harmonic  XAjould  be  apnroximatelp  20  j  too  large),  the  comnarison  should  be 
viewed  as  qualitative.  In  addition,  man^r  of  the  exneriraental  noints  corres- 
nond  to  trace  a.-m'^litudes  less  than  10,3  of  the  total  oscillator7/  amnlitud es , 
and,  therefore,  the  accxiracy  of  these  results  is  uncertain. 


Flight  Tests;  K-3li.i  =  0.13 

Only  rigid-bod.y  first  harmonic  notions  x-xere  used  in  calculating  the  i'I-3ix 
airloads  at  an  advance  ratio  of  O.lO.  These  calculated  airloads  are  compared 
in  Figure  10  with  the  measured  airloads.  There  is  no  measured  eleventh  har¬ 
monic,  so  only  the  comnuted  distribution  is  shox-ra  for  this  harmonic.  Again, 
stall  effects  and  Mach  number  effects  were  neglected.  Agreement  between  the 
comnuted  and  measured  airloads  certainly  leaves  something  to  be  desired,  but 
it  is  again  noted  that  the  first  harmonic  airload  moment  is  not  zero  for  the 
coiBDuted  distribution.  The  calculated  steady  load  is  less  than  the  measured, 
which  is  the  inverse  of  the  situation  found  ■with  the  wind  tunnel  model. 
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Flight  Tests;  H-3iu  “0.29 


For  the  11-3^'  at  an  advance  ratio  of  0,29,  the  calculated  airloads  are 
co’Tirared  vith  measured  airloads  in  Figure  11,  0nl3'’  rigid-body  first  harmonic 

motions  were  included  in  the  comnutations  shown.  Again,  Mach  number  and 
stall  effectvS  were  neglected. 

It  vas  nointed  out  oreviously  that  the  experimental  data  presented  here 
renresents  onlj^  a  ver^r  small  ’portion  of  that  obtained  during  the  IRSa  tests, 
Man;''  conditions  other  than  steadjr  state  were  investigated  during  each  flight. 
This  obviously'-  required  tiie  judicious  setting  of  data  channel  gains  to  insure 
that  data  were  obtained  for  flight  maneuver  conditions,  Unfortunateljr,  these 
practical  limitations  restricted  certain  of  the  blade  motion  and  strain  trace 
deflections  to  relativeljr  small  values  when  steadj^'-state  level  flight  condi¬ 
tions  were  floi'm,  ouperimnosed  on  nossible  reading  inaccuracies  are  the  errors 
of  untaown  magnitude  that  result  from  conversion  of  the  strain  and  motion  data 
to  deflections  and  velocities.  For  these  reasons,  and  because  compensating 
errors  might  be  introduced  imnlicitly  through  the  apnroximations  used  in  the 
theorj.'-,  the  calculations  viere  made  neglecting  the  blade  elastic  motions. 
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CONCLUD ING  UGJ I4HKS 


It  is  believed  that  the  results  obtained  to  date  in  this  continuing 
study  indicate  that  accurate  calculation  of  the  periodic  aerodynamic  loads 
on  rotar;;^  wings  is  feasible.  Permissible  approximations  and  assurriptions  are 
not  yet  clearly  defined.  The  effect  of  wake  distortion  may  be  important, 
but  a  good  approximation  for  the  walce  displacement  time  histories  should  be 
found  and  used  before  a  conclusion  is  reached,  Sven  the  number  of  spanwise 
stations  required  for  computation  purposes  may  be  a  function  of  flight 
condition  due  to  the  rapid  loading  changes  encountered  when  sections  of  the 
blades  stall  or  Mach  number  effects  become  important. 

Treatment  of  the  combined  performance,  aerodjTiamic  loads,  and  blade 
response  problems  is  believed  to  be  necessaiy,  Moviever,  the  possibility  of 
considering  the  blades  to  be  rigid  for  purposes  of  calculating  the  forcing 
functions  and  then  deriving  the  elastic  response  and,  subsequently,  super¬ 
imposing  corresponding  damning  loads  cannot  be  rejected. 

The  above  discussion  is  intended  onl3’’  to  indicate  that  the  results 
obtained  are  encouraging  even  when  the  predicted  and  measured  airloads  are 
compared  on  a  most  stringent  basis.  In  some  cases  a  small  phase  shift  would 
produce  a  redistribution  of  the  sine  and  cosine  components  to  give  better 
agreement.  In  additj.on,  the  calculated  values  have  borne  the  burden  of  any 
errors  in  blade  position  measiu’ement,  strain  measurement,  and  the  blade 
physical  properties,  ’'^ind  tunnel  corrections  were  not  applied  to  the  model 


data  and,  of  course,  any  departure  from  steady-state  conditions  in  the  experi¬ 
ments  were  not  accounted  for  in  the  computations. 

Continuing  effort  on  this  program  should  include  work  on  the  following r 

1)  The  noncirculatory  and  separated- flow  forces. 

2)  An  improved  approximation  to  the  wake  displacement  time 
history  (distortions)  and  the  effects  of  the  wake  in  the 
immediate  vicinity  of  the  blade, 

3)  The  satisfaction  of  the  cho-^'hj'ise  boimdary  conditions, 

ii)  The  computation  of  the  flapping  motions  in  addition  to  the 
airloads . 

The  noncirculatory  forces  which  are  always  present  vfhen  the  flow  field 
is  unsteady  and  the  separated-floW  forces  vjhich  arise  when  an  airfoil  section 
stalls  can  be  quite  appreciable  and  should  be  included  in  any  method  of 
predicting  the  airloads  on  a  rotating  wing.  In  connection  with  the  stall 
phenomenon  and  the  reversed  flow,  there  are  the  attendant  problems  of  wake 
position,  its  point  or  area  of  attachment,  and  its  nature  (i.e,,  the  amount 
and  orientation  of  vorticity  contained). 

The  induced  velocities  are  now  computed  at  the  blades  and  thus  are  in 
the  rotating  coordinate  system.  For  an  improved  definition  of  the  wake 
displacement  time  history,  it  is  necessary  to  have  the  periodic  induced 
velocity  distribution  in  the  nonrotating  coordinate  system  on  or  j\ist  under 
the  rotor  disk.  It  is  possible,  with  some  modifications  to  the  present 
method,  to  compute  such  a  velocity  aistribution. 
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The  attached  wake,  because  of  its  proximity  to  the  blade,  can  have  a 
significant  effect  on  the  airloads.  Thus,  a  more  accurate  representation  of 
this  portion  of  the  wake  is  desirable.  This  could  be  accomplished  by  replac¬ 
ing  the  trailing  and  shed  vortex  segments  representing  the  most  recent  section 
of  the  attached  wake  (i.e.,  the  section  corresponding  to  the  last  azimuth 
step)  by  a  distribution  more  representative  of  a  continuous  wake. 

The  present  method  is  capable  of  satisfying  linear  chordwise  boundary.'" 
conditions,  but  there  is  some  question  as  to  the  importance  of  accounting  for 
a  more  general  chordwise  boundary  condition  when  corni^uting  the  airloads  on 
a  rotating  vjing.  This  may  be  important  and  should  be  investigated  further. 

In  its  present  form  the  method  de-''-eloped  for  computing  the  airload 
distribution  requires  that  the  blade  motion  be  specified.  However,  when  the 
measured  motions  are  used  for  a  given  rotor  and  flight  condition,  knoi-m 
constraints  are  not  necessarily  satisfied.  For  this  reason  and  the  fact 
that  in  practical  applications  (i,e,,  preliminary  rotor  design)  the  motions 
are  not  knoTrm,  it  x^ill  be  necessarj’’  to  compute  the  blade  motions  simultaneously 
with  the  airload  distributions.  It  is  suggested  that  this  capability  be 
incorporated  first  for  the  flapping  motion  only  and  subsequently  for  the 
remaining  motions. 
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